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We consider the creation of polarization entangled light from parametric down-conversion driven 
by an intense pulsed pump inside a cavity. The multi-photon states produced are close approxima- 
tions to singlet states of two very large spins. A criterion is derived to quantify the entanglement of 
such states. We study the dynamics of the system in the presence of losses and other imperfections, 
concluding that the creation of strongly entangled states with photon numbers up to a million seems 
achievable. 
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Entanglement of light has mainly been demonstrated 
at the few-photon level. It is a challenging goal to pro- 
duce entangled states involving large numbers of photons, 
approaching the domain of macroscopic light. Here we 
propose a scheme that is based on the non-linear optical 
effect of parametric down-conversion driven by a strong 
pump pulse, where the interaction length is increased by 
cavities both for the pump and the down-converted light. 
Our work is thus related to experiments on squeezing [1] 
and twin beams [2, 3]. Polarization entanglement be- 
tween the quantum fluctuations around two macroscopic 
polarized beams has recently been created experimen- 
tally [4]. 

Here we aim to create entangled pairs of light pulses 
such that the polarization of each pulse is completely 
undetermined, but the polarizations of the two pulses 
are always anti-correlated. Such a state is the polariza- 
tion equivalent of an approximate singlet state of two 
very large spins. It is thus a dramatic manifestation of 
multi-photon entanglement. Starting from a spontaneous 
process, the proposed setup builds up entangled states 
which have very large photon populations per mode, cor- 
responding to strong stimulated emission, and thus de- 
serves the name of an "entanglement laser". 

The basic principle of stimulated entanglement cre- 
ation was experimentally demonstrated in the few- 
photon regime in Ref. [5]. To analyze whether the cre- 
ation of large photon number entanglement is possible 
in practice, it is essential to understand how imperfec- 
tions in the setup affect the entanglement. This requires 
a quantitative measure for the entanglement. We de- 
rive a simple inseparability criterion that is formulated 
in terms of the total spin J and the total photon num- 
ber N: if (J^)/(iV) is smaller than 1/2, then the state 
is entangled. Using this measure we show that strongly 
entangled states of very high photon numbers can be gen- 
erated in the presence of losses and other imperfections. 

Let us now study our system in more detail. The source 
of entangled light is described by a Hamiltonian [6] 
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where a and b refer to the two conjugate directions along 
which the photon pairs are emitted, as shown in Fig. 1, h 
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FIG. 1: Proposed setup for an "entanglement laser". An 
intense pump pulse propagates back and forth between the 
mirrors Mi and Ah- Whenever it traverses the non-linear 
crystal C it creates polarization entangled photon pairs into 
the modes a and b, which are counter-propagating pulses in- 
side the cavity formed by the mirrors M3 to Me. The cav- 
ities, which have to be interferometrically stable, are care- 
fully adjusted such that the three pulses (pump, a and b) al- 
ways overlap in the crystal. The fact that a and b propagate 
in the same cavity automatically synchronizes the counter- 
propagating modes. The number of photons in a and b in- 
creases exponentially with the number of round-trips. They 
can be switched out of the cavity by electro-optic switches 
Sa and Sb- The polarization of each pulse is then analyzed 
with the help of polarizing beam splitters (PBS) followed by 
photo-diodes that give a signal proportional to the number of 
photons. Taking the difference between the photon numbers 
for the two polarizations behind each PBS corresponds to a 
spin measurement. The axis of spin analysis is changed by 
appropriate wave-plates in front of the PBS. 



and V denote horizontal and vertical polarization, and k 
is a coupling constant whose magnitude depends on the 
nonlinear coefficient of the crystal and on the intensity of 
the pump pulse. The Hamiltonian describes two phase 
coherent twin beam sources, corresponding to the pairs 
of modes ah,b^ and ay,hh- In the absence of losses, it 
produces a state of the form 
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where t = Kt is the effective interaction time and 
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All terms in the expansion in Eq. (3) have the same 

magnitude, such that the observed polarization (the dif- 
ference in the number of horizontal and vertical photons) 
will fluctuate strongly. However, there is a perfect anti- 
correlation between the a and b pulses. The state \ip) 
looks the same if the axis of polarization analysis is ro- 
tated by the same amount for the a and b modes. It is the 
polarization equivalent of a spin singlet state [7], where 
the spin components correspond to the Stokes parame- 
ters of polarization, 

= \{aiah- ala^), = ^{ala+ - ala-) 

Jy = ^{ajai - alar), (4) 

The spin components can thus be expressed as differ- 
ences in photon numbers, where a+- = -^{afi±ay) cor- 
respond to linearly polarized light at ±45°, and ai^r = 
-^{dh ± idv) to left- and righthanded circularly polar- 
ized light. The label A refers to the a modes, cf. Fig. 1. 
Analogous relations express in terms of the b modes. 
The total spin satisfies (J^)^ = {J^f + {J^f + {J^f = 
^ ("^ + l)- Number states of the modes au and a^ are 
eigenstates of and of (J"^)'^. The state \n — k)ah\k)a^, 
has total spin j = n/2 and eigenvalue m = (n — 2fc)/2. 

The states IV-"") of Eq. (3) arc singlet states of the 
total angular momentum operator J = -|- for fixed 
3 A = Jb = njl. As a consequence, (■!/'|J^|'0) = also for 
the state of Eq. (2). Losses and imperfections lead to 
non-zero values for the total angular momentum, corre- 
sponding to non-perfect correlations between the Stokes 
parameters in the a and b pulses. Since the ideal state of 
Eq. (2) is highly entangled, one expects that states in its 
vicinity are still entangled. We now present a convenient 
criterion for entanglement: for separable states 
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where J = J-^ -|- a.nd N = Na + Nb- To prove this, 
consider (J^) for a separable state p = J2iPiPi 



One has 



i i i 

> E^^*K(J^)'>i + ((J^)'>i - 2|(J^),||(J^),|] 

i 

> EP't^^J"")')' + ^(J'')')' - (6) 



where (J^)i = Trpfj^, (J^),; = Trcrf etc. Further- 
more ai = ^((J^)2)i + i - i,A = + 

and we have used the fact [8] that | (J) | < ^ (J^) + j — 5. 
The last line of Eq. (6) can be rewritten as 

i i 

+Qi + A] > + ^ V^^a) + ^^b)), (7) 
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where the last inequality follows from y (J^) + I ^ 5 
i (TV) , which is a direct consequence of the relation = 
^ + 1). Since A'' = Na+Nb, this concludes the proof 
of our criterion. Thus every state that has {3'^)/{N) < i 
is entangled. This is a tight bound. There are separable 
states that reach {J'^)/{N) = |, for example the prod- 
uct state |2j)a^|0)a„|0)b^|2j)f,^, which in spin notation 
corresponds to \jA = j, niA = j) O IJb = j, ms = -j)- 

It should be emphasized that our criterion is sufficient, 
but not necessary. There are entangled states that are 
not approximate singlets. Our criterion is specifically 
designed for the class of states under consideration and 
for polarization observables. It has some similarity to the 
entanglement criterion for spin-squeezed states derived 
in Ref. [9]. The quantities (J^) and (N) are simple to 
calculate, such that the effects of various imperfections 
can be studied with ease. We start by investigating the 
effect of loss. 

Loss in a general mode c corresponds to a transforma- 
tion c y/fj c + y/1 — rj d, where d is an empty mode and 
7] is the transmission coefficient. Let us start by assum- 
ing that the modes a^, and a^, suffer an equal amount 
of loss described by r]A, while the b modes have a trans- 
mission r]B- Using Eq. (4) this leads to the following 
transformations : 
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The state before losses, Eq. (2), has ((J^)^) = ((J^)^) = 
-(J^ • J«), {N%) = (Nl) = (NaNb) and {Na) = 
(Nb) = {N)/2, which leads to the following expression 
for the total angular momentum after losses: 



(J^) - {AnfiiJ^f) + -[va{1 - Va) + Vb{1 - Vb)]{N),{9) 

where At] = r]A — t/b- Remembering that (J"^)^ = 
^(^ -I- 1) one sees that the first term in Eq. (9), 
which depends on A77, is of order (A^)^, while the second 
term is only 0{{N)). If one wants to observe entangle- 
ment for large photon numbers, it is therefore important 
for the losses (including detection efficiencies) in the a 
and b modes to be well balanced. More precisely, Eq. (9) 
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together with our entanglement criterion imphes the con- 
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An equivalent requirement was met 



for (A^) of order 10® in the experiment of Ref. [3] that 
demonstrated the strong photon number correlations of 
pulsed twin beams by direct integrative detection. An 
analogous condition can be derived for a difference in 
losses between different polarization modes. If all modes 
suffer the same amount of loss, described by a trans- 
mission 7], then only the second term in Eq. (9) remains, 
leading to a loss- induced correction to the ratio {J'^)/{N) 
of ^^^j-^, taking into account that the losses also trans- 
form (N) into r]{N). This gives a critical transmission 
value = 1/3, above which entanglement is provable 
by our criterion. The entanglement is thus surprisingly 
robust under balanced losses. 

So far we have considered a situation where first 
the ideal state of Eq. (2) is created, and then it 
is subjected to loss. However, in the cavity setup 
of Fig. 1, which is required to achieve high pho- 
ton numbers, photon creation (in the non-linear crys- 
tal) and loss (in the crystal and all other optical ele- 
ments) happen effectively simultaneously. It is conve- 
nient to transform to a new basis of modes given by 
Ci = -^{ah+bv)..C2 = -^iah-by),C3 = -^{a^ + bh),C4 

-^(a„ — bh). In this basis the Hamiltonian (1) becomes 
that of four independent, but phase-coherent, squeez- 
ers, H=i^ ((4)2 - (4)2 - (4)2 + (4)2 + h.c}j . Intro- 
ducing the quadrature operators Xi = -^{ci + cl),pi = 



= 2 (^l^^l ~ ^2P2 - X3P3 + X4P4) + h.C. 
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Writing down the Heisenberg equations for this Hamilto- 
nian, xi = i[H, xi] etc., one sees that (pf), (xf), (xf) and 
(pi) become squeezed exponentially, while the fluctua- 
tions in the conjugate quadratures, (a^i), (pi)) (Ps)) (^^l) 
grow correspondingly. In the presence of losses, the 
Heisenberg equations have to be replaced by Langevin 
equations of the form 



Xl = K{t)xi - Xxi + fxl{t) 
Pi = -K,{t)pi - Xpi + fpl{t), 



(11) 



and corresponding equations for the other modes. Here 
the time dependence of K{t) = Koe~^'^ takes into account 
the loss of the pump beam while A is the loss rate of the 
down-converted light; fxi(t) and fpi(t) are the quantum 
noise operators associated with the losses [10], satisfying 

{Uit)Uit')) = (/pi(t)/pi(t')) = -i{Uit)fpiit')) = 

\S{t — t'). Here we have assumed that the loss rate A is 
the same for all four down-conversion modes ah,ay,bh,by. 
We will discuss the case of unbalanced loss rates below. 
Eqs. (11) can be integrated explicitly, leading to 

x^{t) = e-fo '=(*')'^*'a;i(0) + /* di'eX' '='*"^'^*"/xi(i'), (12) 

^0 



where k{t) = K{t) - A and / K{t")dt" = ^(e"^*' - e"^*). 

t' 

There is a corresponding expression for pi (t) where the 
sign of K{t) is flipped. 

To understand what these results imply for the polar- 
ization entanglement, one can express the angular mo- 
mentum in terms of the quadratures Xi,pi. One finds 

Jz = ^{XlX2 + P1P2 - X3X4, - PiPi) 
Jx = ^{XiX3+PiP2, + X2X4+P2P4) 

-^y = ^{-^iPi + X4P1 - X2P3 + X3P2). (13) 

Introducing the generic notation p for the quadratures 

that are squeezed (which are pi, .^2, .^3,^4) and x for 
those whose fluctuations grow exponentially (which are 
xi,P2,P3,X4), one sees that all terms in Eq. (13) have the 
generic form x - p. and one finds (J^) — 3{{x'^){p'^) — j). 
The total photon number N = ^ J2i{^i ~^Pi ~ 1)' leading 
to 



(£!) 

(AT) 



(X2)(p2) 



(a;2) + (p2)_i- 



(14) 



Fig. 2 shows the expected time development of the 

mean photon number (A^) and the ratio {J"^) / (N) as de- 
termined from Eqs. (14) and (12) for realistic parameter 
values. The experimentally achievable value for k can be 
estimated by extrapolating existing experimental results 
[5] to higher pump laser intensities. A value ofr = Kt = 1 
for a single pass through a 2mm BBO crystal is realistic 
with weakly focussed pump pulses of a few /iJ, which is 
still below the optical damage threshold. The cavity de- 
sign of Fig. 1 including switching elements will have loss 
rates on the percent level. Fig. 2 shows that very high 
photon numbers can be achieved with just a few round- 
trips. If balanced losses are the only imperfection, then 
the entanglement is very strong even for large photon 
numbers, as long as the "laser" is far above threshold, 
i.e. as long as the rate of creation of entangled photon 
pairs is much larger than the loss rate (k/A ^ 1). Note 
that we are interested in the onset regime, far from sat- 
uration (depletion of the pump). 

The photon number (N) is limited by the requirement 
of observing entanglement in the presence of other im- 
perfections. In particular, Fig. 2 shows the effect of a 
difference in the loss rates between the a and b modes. 
Suppose that the modes a^, and ay have one loss rate 
Xa, while bfi and by have a different one A^. Then the 
quadratures Xi , pi no longer diagonalize the system. For 
example, xi and X2 satisfy the coupled equations 



Xl = K{t)xi - Xxi - ^X2 + fxi {t) 

X2 = -K{t)x2 - XX2 - -Y^^ + /2:2(*)> 



(15) 
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where A = ^(Aa + Ab),AA = A^ - As and ^^.^a 
are the appropriate noise operators. There are analo- 
gous coupled equations for the pairs pi and p2, and 
X4, and ps and p^. These equations are diagonal for a 
new basis of modes ^i,TTi that is related to the Xi,pi by 
a small rotation, which for AA ^ k takes the simple 
form: xi = Ci + (AA/4k)^2, a;2 = -(AA/4k)^i + 6, 2:3 = 
^3 — (AA/4k)^4,X4 = (AA/4k)^3-|-^4, and identical equa- 
tions for the Pi in terms of the tt^. In analogy to the 
case of balanced losses, the quadratures £,i,tt2,'!T3 and 
^4 grow exponentially, while the quadratures tti , ^2 , ^3 
and 7r4 become squeezed. Substituting the above ex- 
pressions for the Xi,pi into Eq. (13) one finds that, due 
to the small rotation between the old and new diagonal 
modes, the Ji contain terms that are quadratic in the 
new large quadratures (^1, 7r2, TTa, ^4). This leads to an 
0{{N)'^) contribution to J^. The dominating correction 

to the ratio {3'^)/{N) is ^^^{N), leading to the con- 
dition — < ,^ for observing entanglement. In the 



regime far above threshold, where A ^ k, this is fairly 
easy to satisfy even for very large photon numbers. 

The effects of other imperfections can be studied in 
similar ways. The most important one is a phase mis- 
match between the two twin beams, i.e. a Hamilto- 
nian H = in^a^bl — e'^'^alb\) + h.c. instead of Eq. (1). 
This can be brought to the ideal form by a transforma- 



tion a. 
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C3, C4 
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^'^/'^b\, which is equiva- 



C4. 



This corresponds to 
cos f a^a — sin |p3. 



lent to C3 

a rotation of the quadratures x^ 

P3 sin ^X3 + cos ^P3, and analogously for 0:4,^4. Sim- 
ilarly to the case of unbalanced losses, this gives a cor- 
rection to the ratio (J^) / (N) whose dominant term is 
j^(p{N), leading to a condition (p < ^"^^^ for observing 

entanglement. This means that strong entanglement of a 
million photons can be observed if (p is of order tt/IOOO. 
This level of precision of optical phases is challenging, but 
conceivable. Strong entanglement for smaller, but still 
considerable, photon numbers is correspondingly easier 
to achieve. 

An amplitude mismatch in the Hamiltonian, H — 
iK{a\bl — /aj,&jj) -l- h.c. with / real, leads to a differ- 
ent degree of squeezing for the modes ci , C2 compared to 
the modes 03,04, but not to a rotation of the quadrature 
amplitudes, such that the effect on (J^)/(A^) does not 
grow with (iV) . 

Another relevant imperfection is a birefringence- 
related mode mismatch, corresponding to a Hamiltonian 
H = iK{alj)l — alb\) + h.c, where the spatio-temporal 
modes d and b of the vertical light differ slightly from 
the modes a and b of the horizontal light. In analogy 
to the case of losses, one can show that a mode mis- 
match that affects the a and b modes in a symmetric 
way leads to a correction to {3^)/(N) that does not grow 
with (iV), which implies that the birefringence-related 
walk-off, while important, does not have to be reduced 
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FIG. 2: Time development of the ratio (J^)/{iV) and of the 
mean photon number (A*'). The units are chosen such that 
t = 1 corresponds to a single pass through the crystal. The 
initial photon creation rate kq = 1, the mean downconverted 
photon loss rate A — 0.03 and the pump loss rate A — 0.01. 
After 8 passes (A'^) reaches the range of millions. The ratio 
(J^) / {N) is shown for three different values of the loss rate 
imbalance AA, namely 0, 0.001 and 0.002. 



by orders of magnitude with respect to experiments on 
the few-photon level. As before, an asymmetry leads to 
an 0{{N)) effect. Note that the other major errors that 
we have discussed, including the phase mismatch, are 
also related to symmetry breaking between the a and b 
modes. In general, geometric symmetry between the a 
and b modes should be implementable to very high accu- 
racy for the setup of Fig. 1. 

In conclusion, the goal of producing strongly entangled 
singlet-like states of very large photon numbers seems re- 
alistic with our proposed system. Besides extending the 
domain where quantum phenomena have been observed, 
such states would also have interesting applications, for 
example in quantum cryptography [7]. We would like to 
thank W. Irvine, A. Lamas-Linares and F. Sciarrino for 
useful comments. C.S. is supported by a Marie Curie fel- 
lowship of the European Union (HPMF-CT-2001-01205). 
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